Dynamic Scaling Phenomena in Growth Processes 



Mehran Kardar 

Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 



Abstract 

Inhomogeneities in a deposition process may lead to formation of rough surfaces. Fluc- 
tuations in the height h(x.,t), of the surface (at location x and time t) can be probed 
directly by scanning microscopy, or indirectly by scattering. Analytical or numerical 
treatments of simple growth models suggest that, quite generally, the height fluctuations 
have a self-similar character; their average correlations exhibiting a dynamic scaling form, 
[/i(x, t) — h(x!, t')] 2 \ = |x — x'l 2 "/ (|x — x'\ z /\t — t'\) . The roughness and dynamic expo- 



nents, a and z, are expected to be universal; depending only on the underlying mechanism 
that generates self-similar scaling. Despite its ubiquitous occurrence in theory and simula- 
tions, experimental confirmation of dynamic scaling has been scarce. In some cases where 
such scaling has been observed, the exponents are different from those expected on the 
basis of analysis or numerics. I shall briefly review the theoretical foundations of dynamic 
scaling, and suggest possible reasons for discrepancies with experimental results. 
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A. Dynamic Scaling 

The growth of films by deposition is clearly of great technological interest. It has 
recently been recognized that such growth processes also pose important issues in the 
physics of non-linear and complex systems. Rather crudely, we can distinguish between 
three types of growth morphologies: 

1. Layer by layer growth is naturally most desired from the technological stand-point. 
From a theoretical perspective, this type of growth is not stable in the presence of inho- 
mogeneities in the deposition beam. However, this is only a statement about asymptotic 
behavior, as it is clearly possible to grow many layers in this mode. 

2. Unstable growth occurs when the selected orientation of the substrate cannot be main- 
tained. The instability is usually manifested by the formation of mounds, or other macro- 
scopic features on the surface. Another type of instability results from growth that is 
controlled by a diffusive field. This instability may lead to formation of fractal aggregates. 
This is probably unrelated to growth by deposition and will not be discussed here. 

3. Self-affine surfaces appear in a growth mode that is intermediate between the above. 
The average orientation of the surface is maintained, but it becomes rough; the amount of 
roughness grows with time and/or scale of observation in a self-similar fashion. 

The most economical way to characterize self-affine roughness is by appealing to a 
dynamic scaling form[l]: If the height of the surface at location x, at time t, is described 
by a function h(x,t), its average correlations may satisfy the scaling form 

g (r, t) = ([h (x, r) - h (x + r, r + t)f) = |if«/ (J^j . (1) 

Similar scaling forms are frequently encountered in the theory of critical phenomena. Fam- 
ily and Vicsekfl] first proposed such scaling in the context of interface growth based on 
numerical results. There is now a firm theoretical basis for such behavior. The dynamic 
exponent z describes the evolution of correlated regions with time: initially different parts 
of the surface are independent; but regions of correlated roughness form over time, their 
size growing as £(£) oc t x l z . In each correlated region, the width of surface grows as the 
observation scale raised to the roughness exponent a. Thus, the overall width of the surface 
initially grows as t 13 with (3 = a/z, until it saturates at L a were L is the sample size. 
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B. Discrete Models 

A large number of numerical models of growth, of varying levels of sophistication, 
have been developed[2]. Here, I shall briefly describe some of the simplest versions. 
Random Deposition: Particles are dropped randomly over deposition sites, and stick to the 
top of the pre-existing column on that site [3]. The height of each column thus performs 
an independent random walk, such that 



where v is the average deposition rate per column, and D denotes its fluctuations. Since 
the columns are independent, the correlation length does not grow with time, £ = 1 and 
1/z = 0. This model leads to unrealistically rough surfaces whose overall width grow with 
the exponent (3=1/2 without saturation. 

Sedimentation: After the dropped particle reaches the top of a column, it rolls down- 
hill to neighboring columns until it reaches a local height minimum, where it stops. The 
smoothening effects of the rolling particles generate correlations in height, and lead to 
more realistic surfaces. Since the motion of rolling particles is more or less diffusive, it is 
likely that the correlation length grows as f (t) oc t 1 / 2 . This value of z = 2 is supported by 
numerical simulations [4], which find in d = 2 (i.e. for a one-dimensional substrate), the 
exponents a = 0.24 ± 0.01 and (3 = 0.48 ± 0.02. 

Ballistic Deposition: Originally introduced for particles moving in a continuous space [5], 
several variants with discretized particle locations have since been developed [1] [2]. The 
deposited particles are now attached at the first point of contact with the growing aggre- 
gate. Since the first contact may be to the side of a high column, this model leaves behind 
a finite density of voids as it grows. 

Restricted Solid on Solid models: There are no voids, and the height is a single valued 
function of position. These models impose a restriction on the maximum height difference 
between neighboring columns [6]. A growth attempt that violates such a restriction is 
rejected. (This can be regarded as a crude form of desorption.) The last two models are 
expected to be in the same universality class, but better numerical values for the exponents 
are found from the latter. The simulations in d = 2 give a « 1/2 and (3 ss 1/3, while for 
d = 3 they result in a m 0.38 and /3 « 0.24. 

"Realistic models" attempt to more closely mimic the physical processes involved in the 
growth of surfaces [7]. These include the detailed motion of the deposited atom, its sub- 
sequent movement on the surface; incorporation on steps or islands, island nucleation, 
and desorption. The main difficulty with simulating all these processes is that they take 
place over a wide range of time scales, making the analysis of large scale roughness quite 
difficult. 



(hi(t)) = vt, 



and 
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C. Continuum Equations 

Rather than describing the detailed microscopic evolution of the surface, continuum 
equations focus on the (hopefully universal) macroscopic aspects of its roughness, e.g. the 
exponents a and (3. The general philosophy is to examine the evolution of coarse-grained 
(hydrodynamic) variables; in the example at hand, the height function /i(x, t). Whereas in 
near equilibrium situations the evolution equation is obtained from variations of an energy 
functional, such an approach is not appropriate far away from equilibrium. Here we follow 
an approach described in references [8] and [9]: The equation of motion is decomposed as 



Particle deposition is described by the first term. Thus rj (x, t) is a random function 
whose mean value gives the average particle flux at x, and its fluctuations represent the 
shot noise in deposition. For simplicity, we shall assume that the noise is uncorrelated 
at different sites and different times, i.e. it is a Gaussian process with 



Surface relaxation subsequent to deposition is described by the functional $. As we 
shall see in the following examples, it can depend on various properties of the height, 
such as its slope V/i, or curvature V 2 /i. We shall assume that the relaxation is local, 
i.e., it can be adequately described by the first few terms of an expansion of $ in h and 
its gradients. Which terms can be included in such an expansion are then determined 
by the underlying symmetries and conservation laws appropriate to the dynamics. The 
basic idea is that any term that is not excluded for fundamental reasons of symmetry 
or conservation, will be generically present. 

We shall now describe continuum equations for the discrete models introduced in the 
previous section. 

Random deposition with no subsequent relaxation corresponds to $ = 0. Integrating 
d t h = ??(x, t) yields 



dh 
~dt 



T](xL,t) +$ [h, V/i, V 2 /i, • • •] . 



(3) 



fa CM)) = v 
(8rj (x, t) 5?] (x', if)) =DS(x- x) 5 (t - if) 



(4) 




(5) 



from which we can immediately obtain 



{h (x,t)) = vt, 



and 



(Sh (x, t) Sh (x', £)) = DtS (x - x') , 



(6) 
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corresponding to j3 = 1/2, and a zero correlation length. 

Sedimentation was originally analyzed by Edwards and Wilkinson (EW) in ref.[10]. They 
concluded that the main relaxational process is proportional to local curvature, leading to 



d t h = ri(x,t) + vW 2 h. 



(7) 



This linear (diffusion) equation is readily solved in Fourier space, leading to 



(|Mk,^| 2 > 



D 



(8) 



UJ 2 + is 2 k 4 ' 



Recasting the above result in real space then leads to the exponents 



a = 



3-d 
2 



(3 = 



3-d 
4 



(9) 



in general space dimension d. In particular, in d = 2, a = 1/2 and (3 = 1/4, while for 
d = 3, the mean-square width grows logarithmically in both space and time. 
Ballistic deposition is generically a non-linear process. There is no a priori reason why the 
relaxation function $ should not depend on the slope V/i. (A direct dependence on h itself 
is ruled out by the translational symmetry h — > h + constant, of the underlying dynamics.) 
By symmetry, the relaxation should be the same for slopes ±V/i, and hence the first term 
in an expansion in powers of slope starts with (V/i) 2 , leading to 



Higher order terms can also be added, but are irrelevant in that they don't change the 
scaling properties. There are several excellent reviews of eq.(10), known as the KPZ 
equationfll]. Since the origin and properties of the nonlinear term are discussed in detail 
in these reviews [12] [13], I shall not elaborate on them any further. In d = 2, the non-linear 
equation can in fact be solved exactly, and leads to a = 1/2 and (3 = 1/3 (as compared to 
(3 = 1/4 for the EW equation), in excellent agreement with the numerical simulations [1] [2]. 
There are no exact solutions in d = 3, but based on simulations we can estimate a 0.38 
and (3 ps 0.24. A characteristic signature of eq.(10) is the exponent identity 



obeyed in almost all simulations. The KPZ equation appears to describe the asymptotic 
behavior of most local, random growth processes. 




(10) 



a + 




(11) 
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D. Conservative "MBE" Models 

We may ask why the KPZ nonlinearity is not present in the model of sedimentation. 
The reason is that it is forbidden by a conservation law. If the growth process does not 
allow the formation of overhangs or voids, and there is also no desorption, then all the 
incoming flux is incorporated into the growing aggregate. This means that the net mass 
of the aggregate, proportional to H(t) = f d d xh(x., t), can only change due to the random 
deposition; the relaxation processes should conserve H(t). This immediately implies that 
the relaxation function must be related to a surface current, i.e. 

$ = -V-j[V/ i ,V 2 /i,---]. (12) 

The KPZ non-linearity is thus ruled out, as it cannot be written as the divergence of an- 
other function. Following some observations of Villain[14], many studies have focused on 
such conservative models in the context of MBE growth. The basic idea is that aggregates 
formed by the MBE process are typically free from holes and defects, and that the des- 
orption of the adsorbed particles is negligible. It is thus argued that, at least over some 
sizeable pre-asymptotic regime, the relaxation processes should be conservative. Some 
examples of such conservative models are discussed in the remainder of this section. 
Surface diffusion currents in equilibrium can be related to variations of a chemical potential 
via j oc — V/U. Since in the equilibrium between two phases, the chemical potential is 
proportional to local curvature (fx oc — V 2 /i), this leads to an equation of motion 

d t h = t/(x, t) - kV 4 /j. (13) 

This equation is again linear, and as for eq.(8) can be solved by Fourier transformation to 
yield, 

The corresponding exponents in real space are, 

z = 4, a = —— i l3 = ^T' ( ) 

This leads to a = 3/2 and (3 = 3/8 in d = 2, and a = 1 and (3 = 1/4 in d = 3. Note that a 
surface remains self-affine, maintaining a well defined orientation, only as long as a < 1. 
The above large values of the exponent a indicate a break-down of validity of the above 
equation. 
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Non-linear "MBE" models have been proposed partly to remedy the break-down of the 
linear equation for a > 1. One such model starts with a non-linear chemical potential 
introduced by Sun, Guo, and Grant [15], resulting in 

d t h = r](x,t)-KV 4 h + ^-V 2 (Vh) 2 . (16) 

Despite its non-linear form, eq.(16) can in fact be analyzed to yield the exact exponents[16] 

7 + d 5-d 5-d 

Z= ^ > Ot — , U = (17 

3 ' 3 7+d v ' 

In particular a = 2/3 and (3 = 1/5 in d = 3. However, in non-equilibrium circumstances, 
there is no reason for the surface current to be derivable from a chemical potential. Remov- 
ing this restriction allows the inclusion of other non-linearities [17], such as V (V/i) , which 
are in fact more relevant. More importantly, non-equilibrium currents should generically 
include a term proportional to V/i, which is the dominant gradient as discussed next. 
Diffusion bias refers to such generic non-equilibrium currents that are proportional to the 
local surface slope[18]. One possible origin of such currents is in the Schwoebel barriers[19]: 
Atoms on a stepped surface are easily incorporated at the step to a higher ledge, but are 
reflected by a barrier towards jumping to a lower ledge[20]. This sets up a net up-hill 
current [14] j = z/V/i, leading to an equation of motion 

d t h = t/(x, t) - v'V 2 h + ■ ■ ■ . (18) 

This equation leads to an unstable growth of fluctuations, and therefore higher order 
terms are necessary to ensure stability. For example, Johnson et al[21] have proposed the 
following non-linear equation 



dh , . „/ u'Vh 
— = rj x,f - V 



7?(x, t) - V ( : _, n2 ) - nV A h. (19) 



The instabilities in this equation develop into a complex array of mounds dubbed as SLUGs 
(Super Large Unstable Growths). Finally, it is also possible for the non-equilibrium cur- 
rents to be oriented down- hill (as in sedimentation), in which case the behavior is the same 
as the EW equation (7) discussed earlier. 



E. Discussion 



An excellent review of the experimental research in this subject is provided by Krim 
and Palasantzas[22]. The observational methods include diffraction (specular or diffuse 
X-rays, RHEED, LEED, HR-LEED, and helium scattering), direct imaging (STM, AFM, 
SEM, TEM), and surface adsorption. A variety of metallic (silver, gold, copper, iron), 
and other (Si, InP, NbN, polymer) surfaces grown under a host of different conditions 
have been examined by such probes. Some of these surfaces exhibit unstable growth, 
while others appear to satisfy self-similarity scaling. However, there is usually no clear- 
cut identification of the exponents with the theoretical models. Some experiments on gold 
and silver give roughness exponents consistent with the KPZ value, but larger values of 
p. Other surfaces give larger values of a, consistent with those of the non-linear "MBE" 
equation (17). The reader is referred to this review article [22] for the details. Perhaps 
the following statements at the end of the review are most revealing of the experimental 
situation: "Over %50 of the experimental work reported on here was published in the 
interval from January 1993 to August 1994. The pace of experimental work is clearly 
accelerating, and rapid advances in the field can be expected." 

Given the discrepancies between experiment and theory, we can also ask if important 
elements have been left out of the analysis. The formalism presented so far deals solely with 
a single coarse-grained variable, the height /i(x, t). Other variables may play an important 
role in the evolution of h. For example, in many cases the roughness is intimately related to 
formation of micro-crystalline grains. Variations in crystallinity have so far been left out of 
the theoretical picture. In principle, one could introduce an additional "order parameter" 
M(x, t) describing the local degree of crystallinity. Surface relaxation may then depend on 
this order parameter, leading to $ (V/i, V 2 /i, M, VM, • • •) . We should then also include an 
additional dynamical equation for the evolution of M. Similar coupled stochastic equations 
have been proposed in the context of surfactant mediated growth [23]. A simple variation 
of ballistic deposition which exhibits complex grain formation as a result of a sublattice 
symmetry breaking has also been proposed[24] in the context of growth of diamond films. 

Equation (3) can be regarded as representing a complex filter, converting the white 
input noise ??(x, t) to the correlated random function /i(x, t) through the action of So 
far, the focus has been on the relaxation function $, assuming that the input noise is un- 
corrected. Not surprisingly, correlations in the input noise lead to more correlated surface 
roughness with larger values of the exponents a and /?[25]. Maybe, in view of the rather 
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large exponents observed experimentally, this point should be further investigated. Start- 
ing with -D(k, uj) = (|r/(k, cu)| 2 ) oc |k|~ 2p , roughness exponents may be tuned continuously 
by changing the parameter p. This is not satisfactory, as it relies on a rather arbitrarily 
chosen exponent for noise correlations. Here I shall propose another possibility which is 
less arbitrary, and may be relevant to some experiments. My choice of noise correlations is 
motivated by the scaling phenomena observed in turbulence, another interesting problem 
in non-linear physics. 

In the case of turbulence, the fluid is stirred at long length scales, setting up a Kol- 
mogorov energy cascade[26] by which energy is transferred to shorter wave-length modes, 
ultimately dissipating at a microscopic scale. In the intermediate (inertial) regime, the en- 
ergy density follows simple power laws. We may similarly assume that the deposition noise 
is correlated over long distances £, such that -D(k) is large only for |k| < l/£. Distances 
intermediate between the atomic scale and I are then analogous to the inertial regime in 
turbulence. Adapting the arguments of Kolmogorov to the KPZ equation leads to 

4 2^1 

^~g, «~3> ^~2' ^ 

in all dimensions. (Because there is no conservation of a corresponding energy for the 
KPZ equation in d ^ 1, the arguments leading to eq.(20) are considerably less compelling 
than those in turbulence. Hence, the above exponents should be regarded as a first ap- 
proximation.) The exponent (3 = 1/2 is similar to that of random deposition. Indeed, an 
early experiment on recrystallization of amorphous GaAs films [27], obtained an exponent 
of (3 ~ 0.50 for highly correlated surfaces. I will thus conclude this brief review with a 
question for the experiments: Is it possible to measure directly the correlations in the 
deposition noise? 
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